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We show that the low-momentum scattering of vortex-antivortex pairs can lead to very long-lived 
oscillon states in 2d Abelian Higgs models. The emergence of oscillons is controlled by the ratio of 
scalar and vector field masses, (3 — (m s /m v ) 2 and can be described as a phase transition in field 
configuration space with critical value j3 c ~ 0.f3(6)±2: only models with (3 < /3 C lead to oscillon-like 
remnants. The critical behavior of the system obeys a power law 0(/3) ~ \j3 — /3c\° , where O is an 
order parameter indicating the presence of oscillons and o = 0.2(2) ± 2 is the critical exponent. 



PACS numbers: 11.27. +d, ll.15.Ex, 98.80Cq 

Introduction. In relativistic and nonrelativistic field the- 
ories, models that exhibit spontaneous symmetry break- 
ing support a plethora of topologically-stable field con- 
figurations or defects [1], static solutions of the nonlin- 
ear equations of motion with properties determined by 
the topology of the vacuum manifold. One of the sim- 
plest examples is the kink in Id A</> 4 models, the result 
of a discrete symmetry breaking: the vacuum has two 
disconnected points and the kink interpolates between 
them. Once formed, kinks cannot be destroyed. Unless, 
of course, they meet an antikink. Low-momentum kink- 
antikink scattering leads to the formation of breathers, 
long-lived time-dependent bound states characterized by 
large-amplitude oscillations of the scalar field [2]. At 
larger velocities, breathers have been shown to form only 
at certain windows or resonances [3]. Although A</> 4 
breathers have not been seen to decay in numerical stud- 
ies [3], they may decay via nonperturbative effects [4], 
radiating their energy to spatial infinity. 

In 2d, the simplest stable topological defect with local- 
ized energy is the Nielsen-Olesen vortex found in Abelian 
Higgs models (AHMs) [1]. These vortices are of great in- 
terest in many areas of physics, e.g., as prototypes for 
studying nonperturbative effects in the Standard Model 
and its extensions, or, in the nonrelativistic limit, as 
Abrikosov vortices in the Landau- Ginzburg theory of su- 
perconductivity [5]. Vortex-like solutions may also have 
a variety of cosmological roles, being formed during a 
phase transition at the GUT scale or thereafter [1, 6] or 
during reheating after inflation [7]. 

Inspired by the existence of long-lived breathers in 
Id Xcj) 4 models, in this Letter we investigate the low- 
momentum scattering of vortices and antivortices in 
AHMs. Vortices and antivortices carry equal and op- 
posite topological charge. This allows for their annihi- 
lation. The only adjustable parameter in the model is 
j3 = A/2<? 2 , the ratio of squared masses of the scalar and 
vector fields. We will show that for (3 smaller than a crit- 
ical value (3 C , vortex-antivortex (henceforth vav) annihi- 



lation does give rise to higher-dimensional breather-like 
configurations known as oscillons [8], also characterized 
by long-lived oscillations of the Higgs field magnitude. 
They were found in 2d [9], 3d [10] and higher-dimensional 
real scalar field models [11]. Their properties have been 
studied in great detail [12]. Long-lived oscillatory solu- 
tions have also been found in the decay of sphalerons in 
the Id AHM [13]. Recently, oscillons were found in 3d 
gauged SU(2) models [14] and gauged SU(2)xU(l) mod- 
els [15] when the Higgs mass is twice the W ± boson mass. 

In the above references, the procedure was to solve the 
relevant equations of motion assuming or not spherical 
symmetry and using an initial profile that approximates 
the oscillon solution. So long as certain conditions 
are satisfied [11], the system relaxes into oscillons as 
it evolves in time. Here, we show that oscillons may 
also emerge dynamically, as remnants of vav scattering. 
(Previous attempts to find oscillons in vortex-vortex 
scattering were unsuccessful [16].) Our finding com- 
plements recent results where oscillons were shown to 
form after a rapid quench in real scalar field models 
with symmetric [17] and asymmetric [18] potentials. 
The fact that oscillons may emerge dynamically raises 
their stakes considerably: they can play a crucial role 
during spontaneous symmetry-breaking in a variety of 
physical systems and situations, from the early universe 
to high energy collisions, and possibly in superfluids and 
superconductors. 
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The Model. The Abelian Higgs Lagrangian density is: 

where T>„ = d^—igA^. The scalar and vector masses are, 
respectively, m s = \f\rj and m v = \f2gr\. Their ratio 
defines the parameter = (m 3 /m v ) 2 — \/2g 2 . In the 
nonrelativistic limit, j3 — 1 defines the boundary between 
Type I (/3 < 1) and Type II (fi > 1) superconductors. 
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the energy is E — j] 2 J M cPxTt, where the Hamiltonian 
density in the temporal gauge Aq = is: 

n= ^E 2 + B 2 )+d t ^d t cj>+V i ^-V i cP+^<f>-l) 2 . (2) 
From cq. 1 , one derives the equations of motion 

V^<j> = (H 2 - if) ; = f = 2 3 Im {tfV v 4>} . 

(3) 

Nielsen-Olesen vortices are solutions satisfying 
4>{r) = e m6 f(r) and A(r)g = — (n/gr)a(r), with bound- 
ary conditions at spatial infinity [r = (x 2 + y 2 ) 1 ^ 2 — ► oo] 
f(r),a(r) — » 1. n is the winding number. At the vortex 
center, /(0) = a(0) = 0. Vortices have a quantized 
magnetic flux &b = 2nn/g. For a vav pair, the net flux 
is zero. Quantities are measured in units of n = 1. 

Lattice Implementation. We work on the temporal gauge 
A = and follow the Hamiltonian implementation for 
lattice gauge theories described in refs. [16, 19]. The 
discrete Hamiltonian is: 

H = J2n^,^)Sx 2 , (4) 

X 

where the density is: 

H =irlir x +Vcl>l-'D<p x + ^E x -E x (5) 

+ ^2 E A* - ^ ~ A Ui + 4) 2 + n4*x), 

and the lattice covariant derivative is T)^ = 
(e~ %g5xA *(t>x+n — 4>x)o~x~ 1 . The Hamiltonian density is 
invariant under that lattice gauge transformations: 

<t> x -> e IX *<j) x ] ir x ^e ix *w x (6) 
^ - A£ + (Xx+ju - Xx)/(gdx); E£^E£. 

The equations of motion are then: 

X 

a: 

for each of the fields 7r Q = {7r, Ei} and 0" = {0,^}. 
The temporal evolution follows a second-order leapfrog 
scheme with St < 5x/^/2. Provided that the initial con- 
ditions satisfy the lattice Gauss's law constraint, 

^ E^-i - £ x) - ifatl (8) 

i 

it will remain satisfied during the evolution. This can 
be implemented by, e.g., setting all temporal derivatives 
to zero as the simulation starts. Energy is conserved to 



0(5t 2 ). We used 5t = 0.02 and 5x = 0.2 and checked that 
the results are consistent for 5x = {0.1 <-> 0.3}. Larger 
values compromise the accuracy of the results. 

Since we use periodic boundary conditions, to prevent 
energy from the initial condition to come back to the 
configuration we use a lattice much larger then the region 
occupied by the vav pair. However, as oscillons are very 
long-lived, we also use an adiabatic absorbing wall [9] 
implemented by adding gauge invariant j(x, y)d t (f> terms 
to the evolution equation, where "f(x,y) has support 
only on a band near the lattice boundary and slowly 
increases from — > 0.1. 

Vortex Antivortex Scattering: Results. We set the initial 
conditions by using an ansatz that approximates a vav 
pair separated by an initial distance So- This is then al- 
lowed to evolve in time under the equations of motion, 
eq. 3, at high viscosity until it settles into the vav con- 
figuration. At large 5q and for j3 < 1 (the case of interest 
to us) we computed numerically the total energy of a 
vortex, finding a good fit for 

E v = 27T/3 1 / 5 , (9) 

which agrees near the region /3 ~ 1 with the result of 
[20]. Also, the contribution in gauge fields to the vortex 
energy, E g = (E 2 + B 2 )/2, can be fitted as 

E g = (10) 

After the initial ansatz settles into a vav pair, we let 
them approach until their centers reach the distance Si = 
Am ~ 1 , large enough to avoid any tail overlap [20]. At this 
point, we set the viscosity to zero and let the pair attract 
and scatter. We then repeat the experiment for different 
values of (3. 

In Fig. 1 we show a sequence of snapshots of an anni- 
hilation process for (3 — 0.08, where an oscillon is formed. 
We plot both the amplitude \(j)\ 2 and the magnetic field 
B z . One can see the outgoing wavefront of radiation 
(center) as the vortices coalesce into the oscillon and its 
oscillating dipole magnetic field (right). 

In Fig. 2 we plot the associated energy (H) 
and Chcrn-Simons (ncs) densities, where ncs — 
(Air)^ 1 e a fj 1 A 1 F al3 . Note that oscillons have a distinc- 
tive signature and thus could have an interesting role in 
theories with current anomalies of the form d^Jg cx FF 
leading to baryogenesis, including during postinflation- 
ary reheating [21], a possibility we are exploring. 

Vortex-antivortex scattering can lead to two possible 
outcomes, depending on the value of (3: either the vav 
pair quickly annihilates, (+) + (— ) — > 0, or it settles into 
an intermediate, long-lived oscillatory state before radi- 
ating its energy to spatial infinity, (+) + (— ) — > (H — ) — > 
0. In order to investigate the process quantitatively, we 
measured the energy in a disk (i?disk) surrounding the 
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FIG. 1: Snapshots of \<f>\ 2 (top) and B z (bottom) for a con- 
figuration with (3 — 0.08. Time proceeds from the initial vav 
pair on the left to the oscillon on the right. 



P = 0.22 10 
A 



I s< E=8.05 

E=3.06— ^ 


P=0.108: 







1000 




200 



400 600 

time 



800 



1000 



FIG. 2: Chern-Simons density (top) and energy density (bot- 
tom) for the same snapshots of Fig. 1 



scattering process. The results are shown in Figs. 3 and 
4. In Fig. 3 we show the energies in a disk of radius 
r = 15 for different values of (3. In Fig. 4 we show the 
total energy and the associated energy in gauge fields for 
j3 = 0.052. The inset shows that there are two character- 
istic frequencies, set by the scalar and gauge field masses, 
m s and m v . For = 0.052, one obtains m v ~ 4.4m s . 

By inspecting Fig. 3, we identify an oscillon whenever 
the energy develops a plateau. To read off a value of the 
plateau energy we examine the energy vs. In t plot. When 
oscillons are present, the curve assumes a "hockey stick" 
shape. We read off the oscillon's energy (E osc ) as the 
value of i?disk at the location of the elbow. [Arrow in inset 
in Fig. 3.] We have observed that oscillon lifetimes are all 
in excess of t osc > 10^ 4-5 ^. For example, for (3 = 0.098, 
T sc > 1.1 x 10 5 . We leave a detailed study of their 
extreme longevity for future work. 

We identify (3 C = 0.13(6) ± 2 as the critical value sep- 
arating the two regimes, as it exhibits the critical slow- 
ing down typical of continuous phase transitions. For 



FIG. 3: Energy within a disk of radius r = 15 surrounding 
the initial vav pair for various values of j3 on a 100 2 lattice. 
The range of j3 is indicated on the top left. There are two 
regimes: for > j3 c ~ 0.136, the pair quickly annihilates, 
radiating its energy to spatial infinity. For /3 < /3 C , the pair 
sheds some of its initial energy and settles into a long-lived 
intermediate state, where almost no energy is radiated. This 
state is the U(l) Abelian Higgs oscillon. Inset: change in 
energy (E) radiated as a function of \nt for = 0.108. Note 
the "hockey-stick" shape. The arrow denotes the oscillon's 
energy, E OBC . 
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FIG. 4: Top curve: total energy for /3 = 0.052 in a disk of 
radius r = 15 surrounding the initial vav pair. Lower curve: 
energy in the gauge fields. Black line is the running average 
of the gauge energy. The inset shows the total energy and 
the energy in gauge fields for a small time interval. The two 
characteristic frequencies satisfy m s ~ \J]3m v . 
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FIG. 5: Phase diagram for the vav — > oscillon transition. Plot- 
ted is the order parameter 0(/3) = E oa c/E v as a function of 
(3. The diamonds are numerical results, while the continuous 
line is the fit of eq. 11. 



j3 > 0.2, well outside the critical region, no oscillons form: 
vav pairs radiate their initial energy to spatial infinity so 
that by t ~ 500, £ disk ^ 0.4 < £ osc . [We neglect here 
some small remnant energy -Edisk < 0.1 — 0.2, that per- 
sists long beyond the oscillon's demise.] 

The quantity E osc /E v naturally lends itself as an effec- 
tive order parameter describing the transition from vav 
pairs to oscillons in field configuration space. In Fig. 5 
we plot E osc /E v vs. (3. The curve is well fitted for 



E osc /E v ~(\f3-p c \T-W. 
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Near criticality the exponent 0.2 will have small correc- 
tions [22] . Given the long integration times needed to ex- 
tract a more precise value for the energy of the oscillons 
around (3 C (akin to critical slowing down, as can be seen 
in Fig. 3), we did not attempt to find these corrections 
here. From the theory of continuous phase transitions, 
the critical point is related to infrared divergences. In 
fact, the effective size of the oscillon increases near j3 c : 
disks of larger radii are needed to capture the full energy 
of the configuration. The associated decrease in the ef- 
fective scalar field mass near the core [~ {A 2 ~ /3) 1 / 2 ] is 
the result of screening by the gauge field. Above criti- 
cality, the mass at the core becomes positive, eventually 
matching the asymptotic vacuum as the vav annihilates. 

Our results show that the annihilation of a vav pair 
in the 2d AHM can be described as a phase transition 
in configuration space, with complete annihilation cor- 
responding to the "symmetric" phase E osc /E v — ► 0, and 
oscillons corresponding to long-lived "broken-symmetric" 
states, with the associated order parameter attaining a 
nonzero value. 



It will be interesting to investigate if oscillons can ex- 
ist in superconductors, how their lifetimes vary with (3, 
and if they can be created in non-Abelian models via, for 
example, low momentum monopole-antimonopole scat- 
tering. 

We thank Noah Graham for stimulating discussions. 
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